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Abstract. We report on measurements of the electrical conductivity in both a 2D triangular lattice of 
metallic beads and in a chain of beads. The voltage/current characteristics are qualitatively similar in both 
experiments. At low applied current, the voltage is found to increase logarithmically in a good agreement 
with a model of widely distributed resistances in series. At high enough current, the voltage saturates due 
to the local welding of microcontacts between beads. The frequency dependence of the saturation voltage 
gives an estimate of the size of these welded microcontacts. The DC value of the saturation voltage (~ 0.4 
V per contact) gives an indirect measure of the number of welded contact carrying the current within 
the 2D lattice. Also, a new measurement technique provides a map of the current paths within the 2D 
lattice of beads. For an isotropic compression of the 2D granular medium, the current paths are localized 
in few discrete linear paths. This quasi-onedimensional nature of the electrical conductivity thus explains 
the similarity between the characteristics in the ID and 2D systems. 

PACS. 45.70.-n Granular systems - 72.80.-r Electrical conductivity of specific materials - 81.05.Rm Porous 
materials; granular materials 



1 Introduction 

Granular materials are ubiquitous in the nature and have 
been extensively studied mostly since the 1990s pQ. How- 
ever, transport phenomena within granular media such as 

■ sound propagation [H[3], heat transfer [H[S] and electrical 
conduction [6] are not so much addressed. 

Electrical conduction within metallic granular media 
is a very old problem. In 1890 E. Branly discovered the 
extreme sensitivity of the conductivity of metal fillings to 
an electromagnetic wave [7]. Although this effect was in- 

[ volved in the first wireless radio transmission near 1900, 
it is still not well understood (see Ref. [8] for a review). 
Except some sporadic works in the 1960-70s |9,10J, the 
electrical conduction in granular media prompts renewed 
interest only since recently. Some works revisited the ac- 
tion of nearby electromagnetic wave (Branly effect) [Til 
[12], others focused on the effect of an electrical source 
directly imposed on the system 8,13,14 . Both perturba- 
tions lead to a transition from an insulating state to a con- 
ductive state of the granular medium. Using a ID model 
granular system, it has been shown that the transition of 
conduction is related to the local welding of the microcon- 
tacts between grains when the applied voltage reaches a 
critical value PIT4]. 
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Electrical conduction within granular media also dis- 
play other astonishing properties: stochastic current fluc- 
tuations [IO l [T5 l fT6 l fT7] . slow relaxations [T7) . percolation 
18, 19j. Some of these electrical properties of metallic gran- 
ular packings can be due to the extreme sensitivity of the 
intergrain electrical resistance to the nature of the grain 
surface (20] . One could even suspect completely irrepro- 
ducible behaviours, fully dependent on the history and 
minute details of the material. Curiously enough, generic 
and reproducible behaviours can be observed in some cir- 
cumstances. 



This paper attempts to describe some of these generic 
behaviours. To wit, two model granular systems will be 
compared: a ID linear chain and a 2D lattice of beads 
under an applied stress. Surprisingly, the voltage/current 
characteristics are qualitatively similar in both experi- 
ments. Here, we also report the first experimental map 
showing the current distribution within a 2D granular 
medium. Note that direct infrared visualizations and nu- 
merical simulations of temperature distribution exist [131 
[2Tj . as well as various 2D models of random resistor net- 
work |22j . This work is a first step towards more complex 
systems, some questions being still open, e.g., the stress 
dependence of the electrical resistance of a granular pack- 
ing [53] . 
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Fig. 1. Sketch of the 2D experimental setup (top view). The 
real size of each side of the hexagonal lattice is 31 beads. 

2 Experimental setup 

The 2D experimental setup is sketched in Fig. [1] Stain- 
less steel beads, 8 mm in diameter (with a ±4/iim toler- 
ance 22); are confined in a hexagonal cell between two 
horizontal plates of Teflon (PTFE) in order to achieve 
low friction between beads and framework. The beads are 
placed on a triangular lattice. The total number of beads 
is 2792, and the size of the lattice side is 31 beads. Three 
sides of the hexagon are fixed, while the others may move 
independently along their normal direction. As shown in 
Fig. [TJ an horizontal stress can be applied either in one 
single direction or isotropically with the help of a feed- 
back loop [3]. The isotropic compression is created by 
compressing in the three directions by means of 3 linear 
motors, whereas the compression in one single direction 
is created with one motor, the two other ones being kept 
fixed. The maximum force applied on each side is limited 
to 80 N. The applied force on each lattice side is mea- 
sured by three static force sensors. In the following, only 
an isotropic stress is applied. Each side of the hexagon is 
electrically insulated, except both faces perpendicular to 
the y-axis (see Fig.[lJ. A DC current source is supplied to 
both electrodes by a source meter (Keithley 2400) which 
also gives a measurement of the voltage drop. The cur- 
rent is injected to the lattice side and not to only one 
bead. During a typical run, a current I is applied in the 
range 10~ 6 < I < 0.1 A, and the voltage U (or the re- 
sistance R = U/I) is simultaneously measured. The cur- 
rent is supplied during a time shorter than 1 s in order 
to avoid possible DC Joule heating. Similar results are 
found in experiments conducted by imposing the voltage 
(10~ 2 < U < 200 V) and measuring the current (or the 
resistance). 

In this 2D experiment, the map of current paths is 
monitored by means of a magnetic field sensor (Barting- 



ton Inst. Mag-03) [25 . Current flow through a conductor 
generates a magnetic field B perpendicular to the axis of 
the conductor ((Ersted experiment) and proportional to 
the current strength (Biot and Savart law). The sensor 
sensitivity is of the order of 500 /LtT. To avoid the interfer- 
ence from the Earth's magnetic field, a high frequency (w 
10 kHz) current is applied to the bead lattice by means of 
the source of a lock-in amplifier. A typical current map is 
obtained as follows. The thick Teflon lid (40 mm in thick- 
ness) is replaced by a thin hexagonal glass plate (8 mm in 
thickness) in order to achieve the magnetic measurement. 
A small stress ~ 1 N is imposed on each side to settle the 
contacts. The glass plate is then tightened on the beads 
network by means of 6 screws and a bar (parallel to the x- 
axis) to avoid a possible buckling of the network at larger 
force. We first apply an isotropic stress on the 2D network 
(< 80 N by side). The system is then allowed to relax for 
one hour in order to avoid possible mechanical drift. A 
voltage is then applied to both electrodes parallel to the 
x-axis (see Fig. [1} , the total current being measured with 
an ammeter (Keithley) . We wait until the current is stabi- 
lized. The magnetic sensor then is placed above the glass 
plate and is parallelly moved along the x-axis (orthogo- 
nally to the current direction) into grooves of the plate, 
from bead to bead along each of the 60 lines of beads. In 
the vicinity of each bead, the value of B is then measured 
with the lock-in amplifier. We have checked that this value 
is proportional to the current averaged within a radius of 
two adjacent beads. One hour is necessary to perform the 
whole current map. 

The one-dimensional experiment consists of a chain of 
stainless steel beads, 8 mm in diameter (with a ±4 /mxi 
tolerance [H]), submitted to a fixed stress. The number 
of beads, N, in the chain can be varied from N = 2 to 40. 
A current, /, is injected to the chain, and the voltage U of 
the whole chain is simultaneously measured (in the same 
way as in the 2D experiment) as well as the voltage, Ui, 
between each beads by means of multimeter/switch sys- 
tem (Keithley 2700 with multiplexer module 7702). The 
ID experimental setup has been already described else- 
where [14] . 

Note that for both experiments, after each cycle of 
current, new contacts between each bead are needed to 
obtain reproducible measurements for the next cycle |14j . 
To wit, the applied force is reduced to zero, and we roll 
the beads which creates new contacts for the next cycle. 

3 Voltage/current characteristics: saturation 
voltage 

The typical electrical characteristic of a ID chain of beads 
is schematically displayed in Fig. [2] which can be sum- 
marized as follows [8IIT4]. For a given applied current /, 
the measured U(I) depends on the history of the imposed 
current. Whenever the current remained below the value 
/ since the last renewal of contacts, the behaviour of the 
contact is designated as "up-characteristic" . On the other 
hand, once a larger current has been injected, the be- 
haviour is "down-characteristic" . Down-characteristics are 
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Fig. 2. (color online) Log-Log schematic view of the U - I 
characteristic for a ID granular medium. The currents II, In, 
and Is define the limits of the different regimes of the up- 
characteristic: linear; nonlinear (see Sect. [7]), logarithmic (see 
Sect.[6| and saturating (see Sect. [3] and [5]) 

always reversible: as long as |J| < I m ax (Imax being the 
maximum current previously applied), the relation U(I) 
is symetric [U(I) = —U(—I)] and remains the same what- 
ever the history of I since the last occurence of I — I ma x ■ 
Depending of the applied force, three typical currents 
II, In, and Is can be defined (see Fig. [2]): 

— If I m ax < Il-i up and down-characteristics are identical 
and linear (see Fig. [2|. 

— If II < Imax < In, up and down-characteristics re- 
main identical but are nonlinear (see Fig. [2]) . Section 
[7] is devoted especially to this regime. 

— For I > Is, the up-characteristic shows a constant 
voltage U S at (see Fig. [2]) proportionnal to the number 
of contacts ~ 0.4 V per contact. The down-characteristic 
is nonlinear, revealing the existence of local welding of 
microcontacts between beads [8,14 : the crowding of 
the current lines within these microcontacts generates 
a thermal gradient in their vicinity when significant 
Joule heat is produced. At high enough current (i.e., 
when U sa t is reached), this process leads to the local 
welding of the microcontacts. The contact resistance 
increases with the temperature of the welded micro- 
contact, which itself only depends on the voltage across 
it. This regime will be discussed in the present section. 

— Finally, for Zjv < Imax < Is , up-characteristic progres- 
sively raises up to the saturation voltage (see Fig. (2). 
Down-characteristics are intermediate between those 
of the two neighbouring ranges. This will be discussed 
in Sec. [H 

The new results about the saturation range concern two 
points: the extension to 2D systems (see below) and the 
extension to AC voltages for the ID system (see Sec. [5J. 

The U - I characteristic of a ID chain of N — 38 beads 
is shown in Fig. [3^ (N being the number of beads between 
the two bead electrodes). Such typical curve has already 
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Fig. 3. (color online) Typical U - I characteristics, (a) ID 
chain of N = 38 beads for an applied force F — 100 N. The 
dashed line U sa t — 15.8 V corresponds to a saturation voltage 
per contact of V a = U sa t/(N + 1) ~ 0.4 V. The solid line 
is the theoretical down trajectory (see Ref. [8lll4j). (b) 2D 
lattice of beads for an applied force F — 50 N. The dashed 
line corresponds to the saturation voltage U sa t — 25.4 V. This 
value gives a measurement of the number of welded contact 
(see text). Different symbols correspond to different maximum 
applied current, Imax- The letters A, B and C correspond to 
the current maps of Figs. |4^, b and c. 



been discussed in Ref. [HUH] ■ Notably, a saturation voltage 
per contact, V = U sa t/(N+1) ~ 0.4 V, is observed due to 
the local welding of microcontacts between beads [EMU]. 
Figure [3]d displays the electrical characteristic for the 2D 
lattice of beads. Qualitatively, it is the same characteristic 
than the ID one. This is quite surprisingly at first sight, 
but this can be understood when looking at the map of the 
current paths in the 2D system: the current appears local- 
ized in few discrete linear paths whatever the strength of 
the injected current (see Fig. [4] and Sec. [4] below). Such a 
concentration of current on chains thus explains the simi- 
larity between the electrical characteristics in the ID and 
2D lattices. Moreover, for the 2D lattice, the 25.4 V satu- 
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ration voltage (see Fig. [3)3) gives a number of 64±1 welded 
contacts when compared to the V = 0.4 ± 0.01V per con- 
tact of the ID chain. This estimate is in good agreement 
with the dimensions of the hexagon (31 beads per side). 
Thus, as in the ID system, the saturation voltage is also a 
measurement of the number of welded contact in the 2D 
system. 

4 Visualization of current paths in a 2D 
granular medium 

Figures. S^,, b and c show the evolution of the current maps 
when the applied voltage increases. These three maps cor- 
respond to the letters A, B and C in the U-I characteristic 
diagram of Fig. [3j3, that is, to the nonlinear and satura- 
tion regime. The current mainly appears concentrated in 
few discrete paths as shown in Fig. |4Jl. When the voltage 
increases, most of these paths are enhanced, whereas few 
ones disappear (see Fig. |4}d). For higher applied voltage, 
the previous paths are strongly enhanced, and a bridge 
between two paths is even created (see Fig. Hb). For the 
highest applied voltage, the total current decreases from 
31 to 20 mA (during the 1 hour measurement) due to 
thermal drift, and the stress decreases of roughly 10%. 
The typical observation of the current paths is generic to 
various stresses within our experimental range, and is not 
reduced to a small range of compression. The effect of 
the stress strength on the current path evolution will be 
described in a future work. 

The current maps of Fig. [4] show that the current 
paths mainly appear concentrated on linear chains. These 
preferred paths are reminiscent of stress concentrations 
(the so-called "force chains") carrying an external load 
imposed on granular systems |26 [I27|I28] . The spatial dis- 
tribution of these force chains depends strongly on the 
type of loading (isotropic or uniaxial compression, or pure 
shear) [28] . Only the isotropic case is described here. Since 
the contact resistances strongly vary with the applied force 
(much more than for an assumed elastic or plastic contact, 
see Fig. 3 of Ref. [H]), this could explain the similarity 
between force chains and "electrical chains" . However, as 
shown below, the contact resistances are also widely dis- 
tributed for a given force. This random distribution of re- 
sistances could also lead to the concentration of current on 
percolating backbones [18J. The chains of maximum cur- 
rent thus would not be those of maximum stress, except 
maybe at very high loading. However, we did not perform 
a simultaneous measurement of force and current maps. 
Since the resistance distribution is broad, we probe more 
the contact network (the least resistance paths) than the 
force network. 

5 Frequency dependence of the saturation 
voltage 

Let us now examine the AC electrical response of the ID 
chain. As far as the saturation voltage is concerned, two 
extreme cases can be considered: 
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Fig. 4. (color online) Visualization of the current path network 
in an 2D hexagonal packing of metallic beads. The applied 
voltage increases from bottom to top (corresponding to a total 
current of / = 3, 11 and 31 mA and denoted by the letters A, 
B, and C in Fig.|3b). An isotropic stress is created by applying 
a 20 N force on each upper side. The size of the lattice side is 31 
beads. The upper part has not been measured (bead number 
in the y-axis > 52). Colorbar scales are different in the 3 plots. 
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— Either the contact temperature can follow the voltage 
variations across a contact. Then the 0.4 V per contact 
corresponds to the peak amplitude of the AC voltage 
across a contact. For a sine like variation in the voltage, 
the rms value of the saturation voltage is OA/^/2 ~ 
0.28 V per contact. The temperature variations occur 
at twice the frequency of the voltage. This doubling 
temperature frequency stems from the dependence of 
the contact temperature on the square of the voltage 
(see below). 

— Or the temperature cannot follow the voltage varia- 
tions due to the specific heat of the material. Then, 
the 0.4 V is the rms saturation voltage across the con- 
tact. 

To discuss the intermediate regime between these two 
extreme cases, let us introduce the influence of the spe- 
cific heat in the thermal budget. We approximate the tem- 
perature T and voltage U fields as spherically symmetric 
around the center of the metallic bridge. Half of the Joule 
heating goes in each bead. The current flowing from the 
left to the right bead, and considering the right bead, we 
define u as the (negative) difference between the potential 
at distance r and the potential at the center. Far from the 
center, T — T a and u = -U/2. At the center (r = 0), 
T = T(0) and u = 0. The total heat flow at distance r is 



\2irr 2 



dT 
dr 



(1) 



where A is the heat conductivity. The total current I is 
independent of r and is 



I = — 



2nr 2 du 
p dr 



(2) 



where p is the electrical resistivity. Joule heating within 
the considered bead at distance less than r is then 



2nr 2 udu 

~ uI = 

p or 

The thermal budget then becomes 



Ar 



, dT r 2 udu 



dr 



p dr 



dr' 



(3) 



(4) 



where C is the specific heat of the bead material. Us- 
ing the relation Xp = LT (see Ref. [HUH]), where L ~ 
2.5 10~ 8 V 2 /K 2 is the Lorentz number, Eq. (0} can be 
written as 



,8T 
Or 



udu T 
Ldr Ar 2 J 



dr' 



(5) 



Integrating Eq. |(SJ) from r = to oo gives 



m -T(or = - U T+ l^l r Cr 



i dT(r>) 
dt 



dr'dr (6) 



or, with T m = T(0), 



C/ 2 
4L 



2T 



C. 



,2 dT{r) 
dt 



dr'dr 



(7) 
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Fig. 5. Comparison between the AC characteristics measured 
at 1 MHz (o) in the chain of = 9 beads, and the model ( — ) 
from Eq. (fTTjl with A = 1.15 um/A. Only the top part of the 
experimental curve is shown here to focus on the saturation 
regime. 



This is the modified Kohlrausch equation due to the sec- 
ond term in the right-hand side of Eq. ([7|) taking the tem- 
poral dependence into account. For details on the Kohl- 
rausch equation at the thermal equilibrium, see Refs. [H 
29 j for generality, or Refs. [HUH] for this problem. The in- 
fluence of this last term is examined in the Appendix A 
for a sinusoidal voltage U = U cos uit. We then find 

Tl -T 2 =^[l + 8{lo) co S (2cjt - 0(c))] (8) 



where 



0(w) = 



a = pI /irU being of the order of the radius of the metal- 
lic bridge (see Appendix A). I = I cos cut is the first (and 
main) harmonics of the current, k and <p being introduced 
in Appendix A. 

As explained in [8l[T4"]. the maximum value of Tj^, is 
fixed to V Q 2 /4L by the softening of the metallic bridge, 
with V — 0.4V. Substituting this variable into Eq. (|8j) 
leads to the rms saturation voltage 



2 

rms 



V 2 

' n 



V, 



y/l + (32a 2 w/25K) 



U 2 

o 

2 ~ 1 + 9(u) ~ ' 1 + v /l + (32a 2 w/25 K ) 

(10) 

Defining 7 rms = I /y/2 and A = 8p/ (5ttV ), the AC current- 
voltage characteristics for one contact then becomes 



U z = V z 

rms o 



1 + + 



(ii) 



Figure [5] shows a comparison between measurements 
at 1MHz (o-symbols) and Eq. (fTTj) (solid line). For this 
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Fig. 6. (color online) Zoom of the U-I characteristics in Fig. [3] 
on a semi-log plot focusing on small currents, (a) ID and (b) 
2D experiment. (— ): fits U oc In/. 

comparison, to = 2irf with / = 1 MHz, and the value 
chosen for A = 1.15 /mi/ A is within 30% of that obtained 
(~ 0.9 /zm/A) from the physical properties of the stainless 
steel beads used in the chain (see Appendix A). V — 
U sa t/(N + 1) ~ 0.42 V is given by the DC saturation 
measurements. The uncertainties in the measurements do 
not allow to claim for a full agreement. However, it is 
clear that the model gives an order-of-magnitude estimate 
for the frequency effects, and a good explanation of the 
progressive character of the saturation in such conditions. 
It also confirms the size of the metallic bridge invoked to 
explain the saturation plateau: namely a =~ 50 nm for 
I = 10 _1 A, and proportional to the current I. 

6 Up-characteristics for I N < I < I s 

Let us now focus on the intermediate range of the DC 
applied current In < I < Is (see Fig. [21). In this case, the 
voltage across the ID chain of beads progressively grows 
up to the saturation voltage. For this current range, the 



characteristic of the Fig. [3^ is shown in a semi-log plot in 
Fig. [SJi. It can be modeled over roughly two decades as 

U(x\nl. (12) 

Only an inhomogeneity of resistances occurs between beads 
along the chain. In the 2D experiment, there are both 
stress [26l[27ll2"8] and resistance inhomogeneities. However, 
a similar logarithmic characteristic is also observed with 
the 2D experiment for this range of current, as shown in 
Fig. [6b. 

We argue that it is simply due to the wide distribu- 
tion of contact resistances between beads. The rough ar- 
gument is the following. When the current is progressively 
raised from zero, the contacts with the largest resistances 
reach the saturation voltage first. Consequently, a further 
increase of the current does not result in voltage increase 
for each of these contacts. On the other hand, the contacts 
with the lowest resistances do not contribute either to the 
total rise in voltage. Thus this rise is due to the contacts 
whose resistance is such that they are close to the satura- 
tion voltage V . This resistance is thus V /I. For a varia- 
tion of current SI, their voltage increase is SU = (V /I)8I. 
By integration, this leads to the log dependence of Eq. (fT2")) 
in qualitative agreement with the observation displayed in 
Fig. El 




Fig. 7. Cumulative distribution function of the resistance of a 
chain of N = 17 beads for 480 realizations (in a semi-log plot): 
before (□) and after (o) the saturation regime (see Fig. [3^,). 
(R) = 38.3 Q before (0.29 Q after), and a R = 154.5 Q before 
(0.22 n after). 



A direct measurement of the cumulative distribution 
of resistances in the chain confirms this broad distribu- 
tion of contact resistances between beads (see Figs. [7] and 
[8]). Figure [7] shows in a semi-log plot the typical cumula- 
tive distribution function of the total resistance of a chain 
of 17 beads during 480 realizations. The distribution is 
broad before the saturation regime, and narrower after 
the saturation. The inset of Fig. [5] shows the values of the 
resistance between 2 beads during 600 realizations before 
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Fig. 8. (color online) Inset: values of the resistance, R, between 
2 beads during 600 realizations before the saturation regime 
(under stress of 200 N). (R) = 28 fi ; a R = 43 fi ; e {ln(R)) = 

13 f2 ( ), and ffj n (m =1.16. Main: Cumulative distribution 

function of ln(R) (centered to the mean and normalized by 
the rms value). Solid line is a Log-normal fit of mean -0.08 and 
standard deviation 1.12. 



the saturation regime. These values are spread on more 
than 3 decades. The cumulative distribution function of 
the logarithms of resistances is well fitted (solid line in 
Fig. [8]) by a log-normal distribution [201130] , Even simpler, 
a flat distribution in resistance logarithms is also a good 
model, which is equivalent to approximate the cumulative 
distribution of Fig. [8] by its inflexion point tangent. 

Such an approximation allows to formalize the above 
argument. Assume that the resistances r are such that 
the probability density of their logarithm, I = ln(r/r ), 
is P (£) = 1/(2^0 for -An <l<l m - This uniform dis- 
tribution P D is normalized by the constant £ rn , and r a is 
a constant with the dimension of a resistance. Then, the 
average voltage per contact is 

U = r Q I / P Q e*d£ + V Q / P a d£ 

J-im JHtZt) 

and 



dU 
dl 



In 



P„r'dt = rJl \— T -e 
. r n I 



As (V /r I) » e- e ™, we find dU/dl = P V /I which 
leads to the expected logarithm characteristic of Eq. [T2J 
Figure [7] also shows that, once the saturation voltage is 
reached, all resistances have the same value V /I m ax (see 
also [HI). 



7 Nonlinear reversible characteristics 

We now examine the part of the characteristics corre- 
sponding to II < I < In ( see Fig- EJ)- While capturing 




U (V) 

Fig. 9. Nonlinear reversible characteristics. The characteris- 
tics are reversible up to the point (U c , Different symbols 
correspond to different applied forces. Arrow means that the 
process is reversible by decreasing or increasing the current. 
ID experiment with TV = 40. 
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Fig. 10. Critical power U c Ic (for which the characteristics of 
Fig. [5] become irreversible) as a function of Ro- Ro is the value 
of the chain resistance at low applied current. ID experiment 
with TV = 40. Same symbols as Fig. [9] 



the main origin of the In/ behaviour for Zjv < I < Is, 
the model of the previous section is oversimplified. It sug- 
gests that the U-I characteristics is linear up to the point 
where the contact with the largest resistance reaches 0.4V. 
This is not true. Figure [9] shows that non-linearities ap- 
pear while the U-I characteristics remains reversible. 

Different applied forces, or even different experiments 
at the same force, yield different values for the irreversibil- 
ity threshold (see Fig. [9|). However, as shown in Fig. IT0| 
the dissipated power U C I C at the irreversibility threshold 
has approximately the same value P c = U C I C for a very 
wide range of initial resistances R a , at least for a fixed 
number of contacts. A similar power-dependent thresh- 
old has recently been reported for 3D samples of Copper 
powder [16, 17 . However, here, we do not observe any time 
evolution of the resistance once the threshold is exceeded, 
as in the case of 3D samples. While the influence of the 



8 



M. Creyssels et al: Some aspects of electrical conduction in granular systems of various dimensions 




Fig. 11. Normalized nonlinear part of the current as a function 
of the voltage U. The logarithmic slope is close to 2.33 ~ 7/3. 
ID experiment with N = 40 beads. Same symbols as Fig. [5] 



dissipated power suggests a thermal origin, we have no 
precise model for this threshold. 

In order to characterize the non-linearity, we plot R I— 
U versus U in Fig. QTJ that is, the nonlinear part of the 
current, normalized by the initial resistance at low applied 
current, Rq. The logarithmic plot stresses a power-law be- 
haviour of this nonlinear part as 



R a I-U ~ U a 



(13) 



The power-law exponent a is always found close to 2.3 ~ 
7/3 whatever the applied force. Let us now use previous 
data from 3D experiments with Copper powder .16], and 
plot R I — U versus U as in Fig. [T3J Here also we find 
a similar power-law behaviour of the nonlinear part with 
roughly the same exponent a than the one found for the 
ID experiment of Fig. [TT1 

The situation for 3D Copper powder samples will be 
clarified in [31] . Let us add a few comments about the 
interplay between the irreversibility threshold occuring 
at a constant dissipated power (P c ), and the saturation 
threshold (U = 0.4V per contact). At the former one, 
U 2 < R P C , which means that, for sufficiently small initial 
resistance R a , it will occur before the saturation. However, 
for large enough R a (small applied forces), it should occur 
after. This is impossible, as our explanation of the satura- 
tion implies the irreversibility. Indeed, in such situations, 
we observed that the voltage U can cross the saturation 
value, and continue to grow until the irreversibility occurs. 
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Fig. 12. Normalized nonlinear part of the current as a function 
of the voltage U. The logarithmic slope is 7/3. 3D experiment 
with Copper powder from Ref. [16]. Symbols corresponds to 
different forces applied to the powder sample ranging from 640 
to 850 N. 



the saturation voltage on the AC frequency confirms the 
model of the metallic bridge, including the order of magni- 
tude of its size. In the 2D system, both the voltage/current 
characteristic and the map of the current paths show the 
quasi-onedimensional nature of the electrical conductiv- 
ity within a 2D granular medium. The knowledge of the 
dependence of the current path distribution on different 
types of loads (e.g., pure shear, isotropic compression, ...) 
will be of primary interest. Another extension of this work 
is to understand the nonlinear reversible behaviour of the 
characteristic observed at very low current in the ID sys- 
tem. The most remarkable feature is that the nonlinear 
part of the current is a power law of the voltage with the 
same exponent that the one found in previous 3D experi- 
ments. 
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8 Conclusion 

We have reported experiments on the electrical conduc- 
tivity in ID and 2D model granular media. The results re- 
ported here serve both to confirm previous interpretations 
and to raise new questions. In both ID and 2D systems at 
low current, the wide distribution of contact resistances 
results in a logarithmic behaviour for the voltage/current 
characteristics. At high enough current, the dependence of 



Our goal is to approximately solve the Kohlrausch modi- 
fied equation ((7} 

with U = U cos ut. 

To wit, we shall use a low frequency approximation 
for T(r) in the integrals. The imposed voltage across the 
contact is U = U cosujt. Thus, u — — (U / 2) f(r) cos cut 
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and / = I Q cos u>t. Substituting these two equations into 
Eq. © gives 

U tt 2 df 

7 ° = r TT" ■ 

p or 

Neglecting the temperature dependence of p, we define 
a = r 2 df/dr, a constant length which is of the order of 
the radius of the metallic bridge. As / = 1 for r ~ oo, 
the above relation integrates as / = 1 — (a/r) which is 
valid for r >> a. For r = 0, / = 0. We shall thus use the 
ansatz: 

/(r) = l--£- = -^- 
r + a r + a 

Let us now limit ourselves to the first harmonics in 2u>, 
that is 

T(r) = ^ [h (r) + 0(u)hi(r) cos(2wt - <j>(u))] (15) 

with h o (0) = 1, 4>{0) = 0, 6>(0) = 1 and 9(oo) = 0. 

There is several approximations here. We assume 4>{uj) 
is independent of r, and h\ independent of w. However, 
4>(lu) is probably larger when r increases, as the volume 
of the involved material (oc r 3 ) is larger. In the same way, 
h\(r) is probably smaller at large r for small tu than for 
large u>. Both effects have the result to effectively limit 
the integral to r smaller than the diffusion length Id- We 
thus take them into account by putting the upper limit 
in r to lr> = y/ k/2u>, with k = A/C, A being the thermal 
conductivity and C the specific heat of the bead material. 
For stainless steel: A = 16.2 Wm^K" 1 , C = 3.9 x 10 6 
Jm- 3 K- 1 , re = 4.1 x 1CT 6 nrV 1 and p = 72 x lfr 8 Qm. 

Within the same approximations, at low frequency, 
Eqs. and ([7]) then become 

rW 2 -r 2 = g(i-/ 2 )(i + cos2^) 

= Tlhl + 2TlhM cos(2 W t) - (T 2 - Tlh\/2) 
Neglecting the last term, we obtain 



ir, 



Keeping only the dominant term in Eq. (|16[) . thus tak- 
ing h a ~ \plajr, that is 



Tl = ^-K = ^Jl~p ■ h x = h /2 . 

Looking only at the 2io terms, Eq. ([Til) then becomes 
9(uj) cos[2u)t-4>{u)] = cos(2u>t) + G(uj)0(uj) sin[2ujt- (f)(u)] 
with 

r e ° h (r) r 

G(lo)=Alu / ^ / Cr^hxir^drdr , 
Jo Ar J 



or 



« Jo r Jo 



It thus gives 



9(u) = 



r' 2 h (r')dr' dr . (16) 



(17) 



r ,2 h (r')dr' 



2J2a 



.5/2 



we finally obtain 



^ 2uj 4a _ 4a /2w 
k 5 D 5 V ft 
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